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We study D-dimensional charged static spherically symmetric black hole solutions in Gauss- 
Bonnet theory coupled to nonlinear electrodynamics defined as arbitrary functions of the field in¬ 
variant and constrained by several physical conditions. These solutions are characterized in terms 
of the mass parameter m, the electromagnetic energy e and the Gauss-Bonnet parameter l We 
find that a general feature of these solutions is that the metric behaves in a different way in D = 5 
and D > 5 space-time dimensions. Moreover, such solutions split into two classes, according to 
whether they are defined everywhere or show branch singularities, depending on (m,£,£„)• We de¬ 
scribe qualitatively the structures comprised within this scenario, which largely extends the results 
obtained in the literature for several particular families of nonlinear electrodynamics. An explicit 
new example, illustrative of our results, is introduced. Finally we allow non-vanishing values of the 
cosmological constant length l\, and study the existence of new structures, in both asymptotically 
Anti-de Sitter and de Sitter spaces. 
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I. INTRODUCTION 

The consideration of extensions of General Relativity 
(GR) containing higher-order powers in the curvature in¬ 
variants is motivated by the fact that they arise in the 
quantization of fields in curved space-time [I] and typ¬ 
ically appear in several approaches to quantum grav¬ 
ity such as those based on string theory [||. Indeed, 
it is expected that the effective Lagrangian of quantum 
gravity resulting from the expansion to the low-energy 
regime will generically include these higher-order curva¬ 
ture terms [3j. In this sense, the finding and characteri¬ 
zation of solutions, in particular, black holes, is of great 
interest to shed light on the understanding of this kind of 
modified, classical gravity theories and hopefully provide 
some insights on the low-energy limit of quantum grav¬ 
ity. However, it turns out that the class of gravitational 
actions made up as functions of curvature invariants such 
as g ilv R} tv , R^R^ and R l xvupR <l,JOL ^ generically give rise 
to fourth-order differential field equations of motion and 
brin g i n ghosts. Nonetheless, as found by Lovelock in the 
70’s [5J (see (&] for reviews on this topic), there is a par¬ 
ticular combination of such invariants in which the field 
equations contain only up to second-order derivatives of 
the metric and the quantization of the linearized theory 
is free of ghost 10- The Lovelock Lagrangian consists 
of a sum of dimensionally extended Euler densities as 

[(D- 1 )/ 2 ] 

Llov = X ( 1 ) 

k =0 

where [z] is the integral part of the number z, D is the 
number of space-time dimensions, c& is the fc-th order 
Lovelock parameter and Lk are given by 
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where is the generalized totally antisymmet¬ 

ric Kronecker delta. For a given D only terms with 
k < {D — l)/2 contribute to the equations of motion, 
while terms with k > (D — l)/2 do not and the case 
k = (D — l)/2 becomes a topological term. The zeroth 
and first order terms in the Lagrangian © corresponds 
to the cosmological constant and the Einstein-Hilbert ac¬ 
tion, respectively, while the second order 

Lgb = R^ P oR^ pG - 4+ R 2 , (3) 

is the Gauss-Bonnet (GB) Lagrangian. This is the sim¬ 
plest nontrivial modification of GR providing modified 
dynamics and second-order field equations as long as 
D > 4. For D > 6 one may consider the next order 
in the curvature invariants in © to provide additional 
modified dynamics, and so on. Finding exact, nontriv¬ 
ial solutions to the field equations of the family of La- 
grangians © is, in general, a hard task, only achieved 
in some particular cases Q. Remarkably, the Lovelock 
family of Lagrangians lies at the crossroad of the met¬ 
ric and Palatini formulations of modified gravity, in the 
sense that the field equations in both formalisms turn out 
to be the same, as opposed to the general case [9j. This 
fact makes Lovelock gravities physically appealing and 
further supports the interest on them. Other approaches 
to higher-order gravity theories with second-order field 
equations, dubbed quasi-topological gravities, have been 
recently considered in the literature jHfl. 

Both in GR, Gauss-Bonnet, and Lovelock gravities 
much interest has been paid to charged, nonrotating, 
black holes. This is due, on the one hand, to the fact 
that in this case exact solutions are more easily acces¬ 
sible and, on the other hand, that classical models have 
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been attracted a great deal of attention in order to model 
the behaviour of charged particles. Historically, the first 
example of such models was the Born-Infeld (BI) La- 
grangian, introduced in the thirties to remove the diver¬ 
gence of the electron’s self-energy in classical electrody¬ 
namics Indeed, Bi-like Lagrangians have been found 
to arise in the low-energy limit of string and D-Brane 
physics [ 0 . which has raised a renewed interest on the 
consideration of nonlinear electrodynamics (NED) in a 
variety of gravitational backgrounds. For instance, elec¬ 
trostatic spherically symmetric black hole solutions of BI 
theory have been studied in the context of GR [Hj, in 
Anti-de-Sitter spaces in several dimensions (l3 |. in GB 
theory 0 and in f(R) models 0 - Other physical mo¬ 
tivations to consider NED models include the effective 
Lagrangian description of quantum electrodynamics ef¬ 
fects 17], the finding of regular black hole solutions in 
GR [IM 0 or the existence of models for which Maxwell 
conformal invariance holds in any dimension [20]. When 
any of these modifications of Maxwell electrodynamics is 
coupled to a particular gravity theory, a different struc¬ 
ture of the corresponding black hole solutions arise, which 
manifests in the structure of their singularities, or in the 
number and type of black hole horizons, both qualita¬ 
tively and quantitatively [2ll - 0 . 

The main result of this paper is to show that, no mat¬ 
ter the physical motivation underlying a particular NED 
model, it is possible to establish general statements on 
the features of electro-static spherically symmetric (ESS) 
solutions in GB theory coupled to NED models without 
explicitly specifying the form of the Lagrangian density 
function , provided that a number of physical reasonable 
conditions on the form of the NED Lagrangian density 
function ( “admissibility”) are assumed. This means that 
NED models are defined here as arbitrary functions of 
the field invariants. These models fall into classes, ac¬ 
cording to the qualitative features of the corresponding 
black hole solutions, the same for all models within a 
given class. This framework largely broadens the class of 
NED models studied in the literature since, as opposed 
to the aforementioned studies, a full classification of the 
corresponding gravitating structures in any such admis¬ 
sible theories coupled to GB gravity takes place in terms 
of a few parameters. No explicit form for the NED func¬ 
tion is needed. More specifically, these parameters are 
the relation between the flat-space ESS energy, e 0. and 
the mass parameter, m, the number of space-time dimen¬ 
sions, D, and the value of the GB coupling constant, Z 2 
(and the cosmological constant length, l\, in asymptot¬ 
ically (Anti-)de Sitter backgrounds). In some cases the 
behaviour of these solutions, when defined everywhere, 
differs largely between the space-time dimensions D > 5 


and D = 5. 

As an illustrative example of our analysis, the third- 
order Lagrangian of quantum electrodynamics (the 
second-order is the well known Euler-Heisenberg La¬ 
grangian of Quantum Electrodynamics fl7l ]') is briefly dis¬ 
cussed. In this work we only consider horizons with the 
usual spherical topology k = 1 . 


II. INTEGRATION OF THE FIELD 
EQUATIONS 

The D-dimensional action for the GB theory coupled 
to NED models reads 



d D Xy/^g (R — 2A) + aLcB) + Sned, (4) 


where k 2 = 87 tGd, with Gd the D-dimensional gravi¬ 
tational constant, A is the cosmological constant, and a 
is the GB constant with dimensions of (length) 2 . As al¬ 
ready mentioned, the GB term is a topological invariant 
in D = 4 and thus it does not contribute to the field 
equations, so in order to obtain modified dynamics as 
compared to GR, in what follows D > 5 is assumed. The 
matter (NED) action is given by 


Sned = J d D x\G^jip(X), (5) 

where ip(X ) is an arbitrary function of the field invari¬ 
ant X = constructed with the field strength 

tensor — d^A^. In D = 4 another field in¬ 

variant can be constructed, namely, — ^F fllJ F* ,lv , where 
p*!iv _ j g ^e ^ ua i t j ie s t r ength ten¬ 

sor. However, due to the dependence on F*^ such an 
invariant cannot be defined for D > 4. In addition, for 
D odd, one might add a Chern-Simons term to the ac¬ 
tion 0, but we shall not consider it here and restrict 
ourselves to the electromagnetic field invariant A'. Exam¬ 
ples of nonlinear actions for the electromagnetic field cou¬ 
pled to gravity and considered in the literature include, 
among many others, Born-Infeld fill . EMI, general¬ 
ized Born-Infeld [0, logarithmic [0, Euler-Heisenberg 
00 , derivative corrections to Maxwell 0], power- 
like and conformally-invariant [0], Hoffman-Infeld [0 . 
Coulomb-like [0, or models leading to regular solutions 
0 . 

For a given (unspecified) p(X) function, and for 
ESS solutions, whose unique non-vanishing component 
is E(r) = F tr , the components of the energy momentum 
tensor 


1 That this energy be finite or not depends on the behaviour of the 
ESS field around the center of the solutions which, in turn, is just 
a consequence of the functional form of the matter Lagrangian 
density, see Sec lIIIl 


T/ = 2 tp x F a pF a v - <*/*>(A), (6) 

(where ipx = are obtained as 
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T\ = T; = 2 yxE 2 - <p{X) (7) 

T°i=Tli =-<p(X),i,j = 2...D-l. 

On the other hand, the line element for a static, spheri¬ 
cally symmetric space-time may be written as 

ds 2 = e v ^dt 2 - e x{r) dr 2 - r 2 dn 2 D _ 2 , (8) 

where d£l 2 D _ 2 = TTy=i sin 2 OjdO 2 is the metric 

on the unit (D — 2) sphere. Due to the source symme¬ 
try T t 4 = T” it can be shown that the (t,t) and (r, r) 
components of the Einstein equations lead to a single in¬ 
dependent metric function that may be written, without 
loss of generality, as g a {r) = e v ^ r ' > = e~ x( - r \ Indeed, since 
the general procedure to obtain the metric function g a (r) 
for a given NED source Tjf in GB theory has been nicely 
described in Refj30| (and employed in several particular 
cases, see 0, [25H27| l. we shall not repeat the deriva¬ 
tion and instead briefly summarize the main steps. From 
the action (j4]) the variational principle leads to the field 
equations 


g , = (D-2)( 0 -3 V (13) 

It is worth mentioning that m takes the same form as 
in absence of gravity. Moreover, it remains also unmod¬ 
ified for the GB theory, as a consequence of the source 
symmetry T* = T”. Such a first integral determines the 
ESS field once the Lagrangian density function is given, 
and this field takes the same form, in the Schwarzschild- 
like coordinate system ©, as in the absence of gravi¬ 
tation in spherical coordinates. The Einstein equations 
in this case can be easily integrated using m, leading 
to (in what follows we redefine h 2 = k 2 /lod- 2 , where 

P-1 

cod -2 = r((g-i)/ 2 j the surface volume of the (D — 2)- 
dimensional unit sphere and, subsequently, drop the tilde 
by notational simplicity) 

777 O ‘r 

9o(r) = 1 - ^3 + {D _ 2 )r -D —3 e-fr, g) + \ (14) 

where m is an integration constant related to the ADM 
mass, M, of the solution as (3l| 


G^+Ag^ — R^v ^ d^vREcxG ^ -}-A g^ u — n (9) 


M = 


(D — 2)lod-2 
167t 


to 


where the correction is given by 


where we have defined 


(15) 


G^f = 2 [RR^ - 2R {ia R a v - 2 R^R^p + 

+ R 0 ^ 1 RvaPy] - -^g^L G B- ( 10 ) 

For the line element dH]) these equations lead to the rela¬ 
tion (5 = (D — 3 )(D — 4)a) 

9 a(r) ~go{r) = ^(1 ~g a (r)) 2 , (11) 

between g a {r) and the solution with a = 0, which corre¬ 
sponds to the metric function that one would obtain for 
the very same problem in GR formulated in D dimen¬ 
sions, i.e. static spherically symmetric solutions of the 
Einstein-Hilbert Lagrangian coupled to NED matter J3 
and with a cosmological constant term. The finding of 
such solution go (r) is obtained by taking advantage of the 
fact that X = E 2 (r, q) does not depend explicitly on the 
metric as a consequence of gttgrr = — 1. Using this, the 
field equations X ^(gjxF' 11 ') = 0, for ESS fields, admit a 
first integral given by 

r D ~ 2 tpxE(r) = q, (12) 

where q is an integration constant, related to the physical 
charge Q of the ESS field as 


pOO 

e ex (r,q)=oo D - 2 R D ~ 2 T°(R 7 q)dR (16) 

J r 

which is physically interpreted as the energy of the 
ESS field outside of the sphere of radius r in absence 
of gravity. In (fT4l) the (Anti-de) Sitter (AdS) radius 
l\ = — < ' D ~ 1 2 a ) ~' 2 ' > parameterizes the cosmological con¬ 
stant term. Once the solution in the Einstein gravity ra 
is known, the above equation m can be easily solved as 



where we have defined l 2 = 2a. In the limit a —> 0 this 
expression reduces to the one of m and thus we re¬ 
cover the solution of the Einstein-NED-A system, while 
in the limit q —» 0 we obtain the solution of Boulware 
and Deser [f|. Note that there are two different branches 
for the solution 113 depending on e = ±1, which comes 
from taking a square root in the resolution of a quadratic 
equation for g a (r). In the limit a —> 0, the “plus” branch 
leads to p+ ~ 1 + - { 2k _ 2) e ex (r ,<?)), which 

possesses an opposite sign for the gravitational mass. It 
has been argued by Boulware and Deser that in the vac¬ 
uum case this branch is intrinsically unstable and the 
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associated graviton becomes a ghost Q, suggesting that 
this branch is physically of less interest (see, however, 
[ 32 I]). On the other hand, the “negative” branch recovers 
the right GR limit, and thus in this work only this branch 
will be considered. 


III. THE MODELS AND THE ESS FIELDS 

Up to now the discussion is valid for any NED model 
and static spherically symmetric electrovacuum solu¬ 
tions. Let us now specify the class of NED models that 
shall be considered throughout this paper. This analysis 
extends the one performed in [33] to the D-dimensional 
case as detailed next. First we restrict ourselves to mod¬ 
els for which the definiteness, derivability and single¬ 
valuedness conditions on the NED function <p(X) hold 
in all the domain of definition of X > 0 covered by their 
ESS solutions (see Ref. 0 ). Such a condition is im¬ 
posed in order to avoid multi-branched Lagrangian den¬ 
sities [35| . which carry potential singularities at the junc¬ 
tion points in the effective geometry as seen by photons 
(see the analysis of [36] in the case of GR). 

The second requirement concerns the fulfillment of the 
positive definiteness of the energy functional in absence of 
gravity which, for a diagonal energy-momentum tensor, 
Tpv = diag(p,p r ,p' 2 - ■ ■ -,p D - 1 ), reads simply p = p r > 0 . 
Note that the weak energy condition (WEC) > 

0 , with a timelike vector, implies in additionp 2 = ■•• = 
Pd- 1, which holds automatically for the NED energy- 
momentum tensor 0 , implying that our models will thus 
satisfy the WEC. Explicitly, this constraint on the energy 
density implies 

P = T* = 2Xp x - <p(X) > 0, (18) 

which must be satisfied everywhere for any field configu¬ 
ration. Though wormhole solutions violating the energy 
conditions (phantom energy models) have been consid¬ 
ered in the literature @ giving rise to a very active field 
of research [37], here we shall not get into such consider¬ 
ations and restrict ourselves to models satisfying the two 
conditions above. 

To extract information from m we first assume that 
the energy density vanishes in vacuum (X = E 2 = 0), 
which implies that <^(0) = 0. By analyzing the inequal¬ 
ity m one easily obtains that three conditions on the 
function <p{X) must hold. From the fact that the ESS 
field can grow to arbitrarily large values, we must have 
tpx > 0(VX 7 ^ 0), implying that <p(X) is a strictly 
monotonically increasing function (excepting at X = 0, 


2 In these models it is usually assumed that violations of the energy 
conditions can occur due to quantum fluctuations, at least at 
some scales. 


where its derivative can vanish). If we assume the La¬ 
grangian density to be defined everywhere, then the con¬ 
dition <p(X) < 0 must hold in the region X < 0. Fi¬ 
nally, in the region X > 0 the inequality (fT51) implies 
that the function ip(X)/-\fX must be positive increas¬ 
ing. If we add to this discussion the first-integral (O it 
follows immediately that these admissibility conditions 
endorse the everywhere definiteness, strictly monotonic 
behaviour, and single-branched character of the (asymp¬ 
totically vanishing) ESS solutions. An important conse¬ 
quence of these admissibility conditions is the fact that 
the function e ex (r, q) in (fTfil) becomes a monotonically de¬ 
creasing and concave function of r, as can be easily seen 
by double derivation of dlfil) and taking into account the 
admissibility conditions established above. 

Let us now study the behaviour of the total energy in 
the ESS field in D dimensions, which is obtained as 

pOO 

e(q)=UJ D - 2 / r D ~' 2 T$(r,q)dr 
Jo 

poo 

= uj D - 2 / dR(2qE - R D ~' 2 ip) (19) 

Jo 

= q^e{q= 1 ), 

where e{q = 1) is the solution of unit charge. The condi¬ 
tions for the finiteness of (12U1) are easily obtained using 
Eq. (fl2ll , and amounts to a simple extension of the results 
derived in Ref. [33], for the case of arbitrary D > 4. 

In order for (l20l) to converge at r —> 00 , we assume that 
the field vanishes asymptotically as (see the definition of 
Tf in Eq.Q) 

E(r) ~ £, (20) 

where /3 is a constant, and we must determine next the 
value of the parameter p > for convergence of the energy. 
Using m the associated behaviour for the Lagrangian 
density function is given by 

<p(X) ~ X a , ( 21 ) 

where a = . Inserting these behaviours in (1201) and 

noting that r D ~ 2 ip ~ E it follows that finiteness of the 
energy implies p > 1. Therefore, one has a < (D — l)/2, 
with a lower bound given by a > 1/2 (corresponding 
to p —^ 00 ), Thus the ESS field approaches its asymp¬ 
totically vanishing values slower (p < D — 2), equal to 
(p = D — 2) or faster (p > D — 2) than the D-dimensional 
Coulombian one. Consequently, a NED model for which 
the energy of the ESS field at r —> 00 in four space-time 
dimensions is finite, leads also to asymptotically finite- 
energy ESS solutions in the D-dimensional case (D >4). 

For the central region r —> 0 there are, for admissible 
models, two field behaviours compatible with the finite- 
energy requirement. In class A1 the ESS field behaves 




5 


as in m (so it diverges at r —>■ 0 ) but now with a > 
(D — l)/2 for convergence of the energy there, which 
implies 0 < p < 1. The Lagrangian density behaves as 
in Eg. (1211) with the same relation between a, p and D. 
This implies that if we have a model whose finite-energy 
ESS solutions belongs to this class in D = 4 dimensions 
(implying a > 3/2), they become of divergent-energy 
for some larger D. An explicit example of this is the 
EH model 17], defined by a Lagrangian density function 
(in D = 4) <p{X, Y) = f + £ (4X 2 + 7 Y 2 ) , £ > 0. It 
can readily checked that this model contains finite-energy 
ESS solutions in D = 4 dimensions, but not for D > 4. 
However higher-dimensional finite-energy ESS solutions 
of this class Al can be obtained by supplementing the 
EH Lagrangian with higher powers of X (see Sec II V Cl) . 

It is worth pointing out that if the Al family is allowed 
to cover the case with p > 1 , this leads to another family 
of ESS solutions (AO class), now with divergent energy. 
The D dimensional Coulomb field (p = D — 2) belongs 
to this family. Given the monotonically decreasing and 
concave character of e ex (r,q) for admissible models, the 
behaviour of all the models in this family corresponds 
to taking e —>• oo in the discussion of Sec|TV] below and, 
as a consequence, its behaviour will be similar to the 
Reissner-Nordstrom-GB one (first studied in the first of 
Refs. [Hi). 

On the other hand, in the class A2 the ESS field attains 
a finite value at the center, behaving there as 

E{r) ~ a — br a ,a > 0, (22) 

while the behaviour of the Lagrangian density becomes, 
by using G2D , in 

<p(X) ~ (yx- a y + A, (23) 

where A is an integration constant and the restriction 
7 = a ~ D + 2 < 1 must hold for finiteness of the energy. 
In this family there are two different behaviours for the 
Lagrangian density functions. If 0 < 7 < 1 (a > D — 2) 
the Lagrangian density takes a finite value around r —> 0 
given by the value A = </?(a 2 ), while for 7 < 0 (a < D— 2) 
the Lagrangian diverges as X —> a 2 . The case 7 = 0 
(a = D — 2) is singular, behaving as 

<p(X) ~ - In (a - Vx) , (24) 

but shows the same behaviour as the 7 < 0 one. The 
Born-Infeld Lagrangian 0 is a well known member of 
this family, with er = 4. 

IV. GAUSS-BONNET-NED BLACK HOLES 

Having discussed the behaviours of the ESS field com¬ 
patible with the admissibility requirement and the finite¬ 
ness of the energy, as well as the consequences for the 


Models 

Field 

Lagrangian 

Finiteness of e 

r —y 00 

E(r) ~ /3/r p 

v(X) ~ x a 

a < (D — l)/2 

Al 

E(r) ~ /3/r p 

9{X) - X« 

a > (D — l)/2 

A2 

E(r ) ~ a — br a 

<p(X) ~ (a - Vxy 

7 < 1 


Table I. The NED models and their ESS solutions, together 
with the condition for finiteness of e in Ea. (l20l) . In this table 
a = (p + D — 2)/(2p) and 7 = [a — D + 2 )/a. 


qualitative behaviour of the function e ex (r , q) in (fT7>l) . we 
proceed to study their black hole solutions within GB 
theory. In this section we shall deal with asymptotically 
flat solutions so we set the cosmological constant term to 
zero (A = 0). As we shall see, all the cases with D > 5 
posses the same structure of horizons, but the behaviour 
of the metric in some of the cases for D = 5 shows dif¬ 
ferences, leading to some new structures. This is due to 
the fact that in the equation of the horizons (if any), i.e. 
the solutions of gtt = 0 , for the different configurations 

jf^Zex(r h ,q)=m- l -^r°- b (25) 

the right-hand-side at m —> 0 takes the value m for D > 5 
but m — 1^/2 in the D = 5 case. Now, by taking into 
account the monotonic and concave character of e ex {r 1 q) 
for admissible models, the horizons can be obtained from 
the cut points between the curve q) and the 

beam of curves m — r D ~ 3 — (^ 2 / 2 )r £> “ 5 , corresponding 
to different values of m once / 2 is fixed (see Fig[l]). This 
leads to an immediate classification of the number of pos¬ 
sible horizons. Indeed, since the large-r^ behaviour of 
this beam is governed by the term — r^~ 3 , this implies 
that it can be at most two cut points between the curves 
in Eci. (E5l) . regardless of the sign and value of Z 2 . How¬ 
ever, a new feature arises here as compared with the GR 
case, since the term 

, x 2 ll ( 2 k 2 \ 

a ( r ) = 1 +y m ~ J^2 £ ex{r,q)j , (26) 

inside the square root appearing in 03 can become neg¬ 
ative, leading to a complex metric for radius smaller than 
a certain rg, which is the solution of the equation 

r-r 1 + 2 ll 0 (m - -A__ £CT (rs, q)J = 0, (27) 

where l 2 o is the value for which a'(rs ) = 0 is satisfied, 
corresponding to 

(D -1)(D- 2) + 4 K 2 l 2 a Tl{rs, q) = 0, (28) 

This is a new kind of singularity, dubbed branch singular¬ 
ity (BS), arising at a non-vanishing horizon radius. This 
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fact splits the solutions into those defined everywhere 
and those exhibiting branch singularities, according to 
the existence or not of a value Z 2 0 . On the other hand, 
the asymptotic behaviour of the metric compatible with 
the finiteness of the energy there (see Eqs. m an d (GB) 
is given by 


9a(r) 


1+ al 1 + £ 


1 + 


2 li 


,D — 


r((' 


4k 2 (D — 2)/3q \]V2 


(p — l)(p + D — 2 )r p 


r)n 


(29) 


so for asymptotically coulombian fields (p = D — 2 in 
Eg. (BUI) ') we recover the GB-Reissner-Nordstrom solution 
obtained by Wiltshire in Ref. jHI . 

For all the GB-NED solutions the sign of the quantity 


9k 2 

Z(q)=m- w -^e(q), (30) 

as well as the one of Z 2 , becomes essential for their proper 
characterization. Let us now characterize the Gauss- 
Bonnet-NED black holes by analyzing first those whose 
metric is defined everywhere. 


A. Solutions defined everywhere 

1 . Z(q)>0,ll>0 

As seen in FigJT] the horizons (if any) of the differ¬ 
ent configurations are obtained though the cut points 
between the beam of curves m — — (Z^/2)rf ,_5 and 

the (monotonically decreasing) curve jpz 2 £ex ( r i q)- From 
the positivity of Z 2 any curve of the beam in the right- 
hand-side of Eq. (E5l) for D > 5 is a strictly monotonically 
decreasing function, starting from to > £ ex{i", q), and 

thus there is a single cut point with e ex (r, q) in this case. 
For D > 5 the metric at the center behaves as 


g a (r ~ 0) 1 - 


(2Z 2 S(q)) 1/2 

Z 2 r^ 


(31) 


and, consequently, it diverges there to — oo, which, to¬ 
gether with the asymptotic behaviour in (129[l . confirms 
that there is a black hole solution with a single horizon 
rh, given by the solution of the equation (B51) . 

However, for D = 5, at r —> 0 the above beam of 

curves takes the value ^f-e(q). When > holds, 

there is a single cut point (see Fig|T|) between the curve 
-ppz 2 £ ex{r, q) and the beam of curves in the right-hand- 
side of Ea. BSl) . while two, one (degenerate) or none ap¬ 
pear when < holds. The metric at the center takes a 
(finite and < 1) value given by 


g a (r —f 0) ~ 1 — 


(2££(g)) 1/2 

ll 


(32) 


The derivative of the metric there depends on the central- 
field behaviour of the ESS fields. For the Al class the 
leading behaviour term is given by (0 < p < 1) 


/ / /~\\ yet 

9a[r ^ ~ ~ { P + 3)(2Z 2 E(q , )) 1 / 2 r p 


(33) 


(where g' = dg/dr ) which diverges to —oo. In contrast, 
for the class A2 this derivative takes a finite (negative) 
value given by 


9a( r -l 0) - 


4 n 2 qb 

3(2Z2E(g))i/2+-'-- 


(34) 


Consistently with the previous considerations, if E(g) > 
Z 2 /2 then the metric around the center is negative and 
there is a black hole with a single horizon. Extreme black 
holes may arise (see Fig[l]) in the case 0 < E(q) < Z 2 /2 
(now the metric at the center is negative). They are 
defined by the conditions g a {fh ) = g' a ( r h) = 0, which 
lead to 


g r hextr( < l)^ 1 o( r hextriQ) — 1- 

Their mass is given by 


(35) 


Tllhextr(q) — ^ \ r hextr + Z Q + 2k q&(rhextr , <?)] , (36) 

where an integration by parts of the term s ex (r,q) 
has been performed, and the quantity ^{rhextr, q) — 
f rh t E(r,q)dr is the electric potential measured at in¬ 
finity with respect to the horizon. For any ESS field of 
the form (1201) the potential at infinity vanishes and then 
we have 4> = A 0 (r hextr ,q). 

In this D = 5 case, besides extreme black holes (m = 
iTi ex tr{q ))j there may be black holes with inner and outer 
horizons (to > m ex tr(q )), or naked singularities otherwise 
(to < ni ex t r (q))- Let us mention that for the critical 
value of the mass m = ^-e(q), which corresponds to 

3 

a critical charge q crit = ( 2 k '*1(Z=i ) ) 4 (so E = 0, see 
Eqs. m and (l30l) l. the metric at the center is g a (r = 
0) = 1. Such solutions are on the verge of becoming into 
those of subcases with branch singularities (IIV B 21) and 
(1IV B 3D (see subsection IIV Bl) when to is decreased. 


2. E < 0; ll o < ll < 0 

In this case the beam of curves in Ea. (B5l) can cut the 
curve e ex (r,q) two, one (degenerate) or zero times, re¬ 
gardless of the value of D , since to — Z 2 /2 — r 2 | r _K) is 
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y 



Figure 1. Qualitative procedure to obtain the number of hori¬ 
zons in admissible models and solutions defined everywhere 
(see Sec lIV A)) . The solid curves represent y = £ex(r,q) 
for case A2 (finite slope at r —» 0) and A1 (divergent deriva¬ 
tive). Each curve begins at the value jyz 2 £ (q) (these curves 
have been normalized in this figure to some value s(q), by 
simplicity), see Eg. (12511 . The cut points with the beam of 
(dashed) curves m — 1^/2 — r 2 (with 1^ > 0, and different 
values of m) give the horizons of the different configurations. 
Each curve of this beam takes the value m when D > 5 and 
m — la/ 2 for D = 5. The case with l&< 0 would show a posi¬ 
tive slope at r — > 0 for the beam of curves, but the number of 
cut points with y are the same as in the > 0 case. For BS 
solutions only cut points corresponding to r > rs are relevant 
for the existence of horizons. 


gence of some of the curvature invariants such as 
R, R^R^, R^uapR^ va ^■ This is an unavoidable feature 
for the models studied here (classes AO, Al and A2), 
resulting from the admissibility requirement. However, 
some models studied in the literature in the case of GR 
[n| have been found to support singularity-free gravitat¬ 
ing solutions. In such a case the ESS field at the cen¬ 
ter vanishes and, as a consequence of Eq. (fl2l) . neither 
the ESS held nor e ex (r, q) are no longer monotonically 
decreasing and the associated Lagrangian densities cor¬ 
respond to multivalued functions [35|. As a consequence 
the energy density may become non-positive definite. Let 
us also note that another Lagrangian with vanishing (at 
the center) ESS solutions is the Hoffman-Infeld model 
studied by Aiello et al [26], given by a Lagrangian den¬ 
sity of the form 

Lhi = 2/3 2 (l — ri(X) — log(X)), (39) 

where y(X) = ^|j(l — y/1 — X//3 2 ) -1 . These solutions 
in GB gravity present a double-peak behaviour for the 
temperature function. Such classes of ESS solutions vio¬ 
late the admissibility conditions of section Hill and shall 
not be considered here. 


B. Solutions with branch singularities 


always below £(<?)■ For this class of solutions an ex¬ 
treme black hole is formed for the tangent cut point to 
the beam (1551) . leading to the extreme black hole radius 
r ex tr(q), obtained as 


K r hextr'^'a i r hextn q) 


(D-2)(D-3) , 

2 ^hextr 

[D — 2)(D — 5) ;2 

A 


(37) 


of which Eg. (1551) is a particular case. As in subcase 
(IIV A II) with D = 5, the associated extreme black hole 
mass in these solutions for any D can be obtained by 
replacing the condition (1571) into Eg. (1551) . leading to 


lTlhextr(q) ~ JJ _ ^ 1 /hextr T la r hextr + <1 ' (t hextr, <l)\ ■ 

(38) 

However, although the structure of horizons is the same 
in the D > 5 and D = 5 cases (naked singularities, 
extreme black holes or two-horizon black holes) the be¬ 
haviour of the metric at the center is not: as r —> 0 the 
metric is obtained as in Eo. (15Tl) . so it diverges there to 
Too when D > 5, while in the case D = 5 it takes a finite 
value given by Eo. (1551) which, as opposed to the previous 
subcase, is now larger than one. 

All these everywhere defined solutions have a cur¬ 
vature singularity at the center, due to the diver- 


As already discussed, in some cases the equation of 
the horizons (1551) is defined only beyond the singularity 
radius rs and thus in our procedure of obtaining the hori¬ 
zons the cut points between the curve y = q) 

and the beam of curves in (1251) makes sense only for 
r > rs- There are three cases to be analyzed separately: 


1 . £ > 0 ; < 0 

At the BS point given by the Eg. (1571) the metric in this 
case takes the value 


g (r s ) = 1 + vf- < 1, (40) 

while the leading terms of its derivative there are given 
by 


g' a {r S ) ^ - 


T,'[rs)r s 


2 (D - 2)11 (l + f “r(m - ^ ex (r s , q))) 


1/2 


(41) 

(see Eo. (l551) for the definition of a(rs))- In the present 
case the derivative in m is positive as we approach 
r — rs —> 0 + since a'(rs) < 0. Now the number of cut 
points between q) and the beam of curves in 
















(l25l) depends on the ratio r|/Z^, ^ — 1. Indeed, when 
0 > r s /l\ > —1 the metric is everywhere positive and 
thus we are led to a naked singularity, while if r| /l^ < — 1 
we have a black hole with a single horizon. The location 
of such a horizon merges with the one of the BS point in 
the limit —» —1 . Note, however, that the condi¬ 

tions determining the existence of BS (the values of l ao 
and £(<?)) are not independent; indeed they are related 
by Eg. (1251) . since the NED energy-momentum tensor en¬ 
ters into both of them. 


Parameters 

Range 

Horizons D > 5 

Horizons D = 5 

E(«) >0,£ >0 

DE 

1 

2 , 1(e), 0, or 1 

M 

A 

o 

jf^: 

o 

A 

if£ 

A 

o 

DE 

2 , 1(e), 0 

2 , 1(e), 0 

E > 0; <0 

BS 

1 , 0 

1 , 0 

E < 0;Z^ > 0 

BS 

2 , 1(e), 0 

2 , 1(e), 0 

M 

A 

o 

if£ 

A 

sf^: 

o 

A 

o 

BS 

1 , 0 

1 , 0 


Table II. The GB black holes for admissible models, and the 
corresponding number of horizons. In this table we use the 
labels: DE: defined everywhere, BS: branch singularities, e: 
extreme 


2. E < 0; ll > 0 

As in subcase (IIVB ID the metric at the BS point rs is 
finite, but now g(rs) = 1 + r|/Z^ > 1 and its slope there, 
given by Eq. m, becomes negative in this case. Conse¬ 
quently the associated gravitational configurations cor¬ 
respond to black holes with two horizons, extreme black 
holes or naked singularities, depending on the value of 
the mass as compared with the extreme one, m ex t r given 
by Ea. ([38]) . 

3. E < 0; ll < ll 0 < 0 

In this case the BS condition (1271) is satisfied twice, i.e., 
rs < rs>, and, consequently, the metric is only defined 
for r < rs and r > rg/. Since e ex (r, q ) is a monotonically 
decreasing function, and larger than m at r —> 0, as we 
increase the radius, the function a(r) becomes negative 
in an interval (when < Z^ 0 ), before becoming positive 
again. The metric diverges at the center to +oo as a con¬ 
sequence of Eq. m- As pointed out in the first Ref. [l5| , 
for an observer in the asymptotic region of space-time 
only the region beyond the outer singularity radius is 
physically accessible, which implies that in this case we 
are led to single-horizon black holes or naked singulari¬ 
ties, depending on g a {rs') ^ 0. 

Consequently we see that black hole solutions in GB 
theory coupled to NED models supporting finite-energy 
ESS fields [see Table an somewhat interpolate between 
GB-Reissner-Nordstrom-type solutions (e(q) oo) with 
a timelike singularity Csubcases lIV A 21 HV B 2l and lIVB 31) 
and GB-Schwarzschild-type solutions (e(q) = 0) with 
spacelike singularities (subcases IIV A II and IIVB II) . 


C. A particular model 

As an example of the above solutions, which captures 
the main features of the analysis performed here, let us 
consider a model belonging to the Al class of solutions, 
given by the Lagrangian density 


with a\ and a .2 being positive constants, in order to sat¬ 
isfy the admissibility conditions. In D = 4 the two first 
terms in (l42l) define the Euler-Heisenberg effective La¬ 
grangian (in this case an additional term in Y 2 must be 
added to (142)) , which vanishes for ESS solutions in D = 4 
and cannot be defined in D > 4) of quantum electrody¬ 
namics (QED) B3- The EH model has been shown to 
contain finite-energy ESS solutions in D = 4 [34|. How¬ 
ever, for D > 4 the term in X 3 must be added in order 
to keep the energy finite. For D > 6 higher order terms 
in A' must be added to flU for this finite character of 
the energy to hold. Let us note that these terms con¬ 
taining higher powers in the field invariant X arise from 
a low-energy expansion of QED, once the heavy degrees 
of freedom are integrated out in the path integral of the 
original action [38] (see also [39|l. By simplicity, here we 
shall restrict ourselves to the cases D = 5,6, which will 
allow us to illustrate the different structures found here. 

Using in the field behaviour at r —> oo and as r ~ 0 
becomes 

E(r^oo,g)~^; E(r~0,q)~ -ir, 

( 43 ) 

and thus we are dealing with an asymptotically coulom- 
bian ESS field, belonging to the class Al as r ~ 0 if 
D < 7, as can be seen from (1751) and the considerations 
of section III. The behaviour of the associated energy 
density in those limits is given by 

r 0 °(r-»oo,g)~ 2r2 f 0 _ 2) 

6 

T o °(r~0,g)~5(^] (44) 

\3a 2 / r s 

The total energy can be calculated by integrating by 
parts in (l20l) and using y = E(r) as the integration vari¬ 
able, leading to (an additional term, vanishing for the 
finite-energy ESS solutions, has been omitted in this for¬ 
mula) 


<p(X) = y + a\X 2 + a 2 A 3 , 


(42) 


£(<?) = 


A(jJD-2 

D- 1 


r y{oo) 


dy 


ly(0) (y ■ <px(X = y 2 )) D - 2 


(45) 
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gct(r) 



Figure 2. Behaviour of the metric function for the model 
(1-421) with the parameters ai = a.i = q — 1, Z 2 = 2, in the case 
E(g) > 0, > 0. In D = 5 the energy is e(q = 1) ~ 2.508 and 

there are four structures: (I) m = 2.52: naked singularity, (II) 
m ~ 2.841: extreme black hole, (III) m = 3.45: two-horizons 
black hole and (IV) m = 4.2: single-horizon black hole. For 
D = 6 (dashed curve, m = 4.25) we find single-horizon black 
hole. All solutions are asymptotically flat. 


which, for the family (Hi?l) reads 

4wd_ 2 o-i [°° dy 

e(q) = ~ F :—- — > 

U 1 “' 0 {y{\ + 2 aiy 2 +3a 2 y 4 )) Tr ^ 

(46) 

whose value can be obtained once the model parameters 
ai,a 2 are given. It can be easily checked that this en¬ 
ergy is always finite for D < 7, regardless of the values 
of the constants 01 , 0 : 2 , as expected. For this model we 
failed to explicitly work e ex {r,q) out but, however, it is 
possible to characterize numerically the different grav¬ 
itational configurations for this mode@, using Eq. ra¬ 
in Fig l2] we have plotted the behaviour of g a (r) for the 
model parameters a\ = 0:2 = 1, unit charge (for which 
e{D = 5,q = 1) ~ 2.508 and e(D = 6,q = 1) ~ 3.446), 
la = 2 and £ > 0. As expected, the structures for this 
case are those discussed in subsection llV AI we see that in 
the D = 5 case the metric around the center is finite and 
its slope diverges there to — 00 , leading to four classes of 
configurations (naked singularities, extreme black holes 
and black holes with two or a single non-degenerate hori¬ 
zon), while for D = 6 the metric diverges to — 00 as r ~ 0 
and there is a black hole with a single horizon. 

In Fig[3]we have plotted the metric behaviour for the 
second class of solutions defined everywhere (see subsec¬ 
tion |TVA]), corresponding to £(g) < 0 and < Z 2 = 
—0.05 < 0, with model parameters a\ = 02 = 1 and a 
unit charge. For D = 6 the metric at the center diverges 
to +00 and one finds naked singularities, extreme black 


3 By simplicity, in this example we take units in which wd— 2 ^ 2 = 
1. 


g a (r) 



Figure 3. Behaviour of the metric function for the model (1421) 
with the parameters ai = 02 = 1 and q = 1 in the case E(q) < 
0 ,£ 0 < la < 0,1a = —0.05 and for D = 6 . As in Fig[5]the 
set of solid lines (diverging to +00 as r —» 0 ) corresponds to 
different values of m, leading to naked singularities, extreme 
black holes or two-horizon black holes. The plot in the small 
frame shows the behaviour for D = 5. In this case the metric 
as r —> 0 is finite but the associated structures are the same 
as for D = 6 . 


holes or two-horizon black holes. Similar structures are 
found for D = 5 (see the small plot in FigO but in this 
case the metric at the center is finite and larger than one. 
The BS solutions of subsection llV Bl for this model can be 
obtained in a similar way. All these results for this par¬ 
ticular model are in complete agreement with the general 
statements of the previous sections. 


V. NON-VANISHING COSMOLOGICAL 
CONSTANT A ^ 0 

Let us consider now the effect of a non-vanishing value 
of A in the metric function (ED- The interest on black 
holes in AdS spaces is mainly motivated by the AdS/CFT 
correspondence [il}, which establishes a relation between 
the thermodynamics of a black hole and the conformal 
dual field theory lying on the boundary of the AdS space. 
It should be noted that these nonlinear corrections to 
the dynamics of the electromagnetic field are expected to 
modify the physics on the CFT side, which could have an 
impact, for instance, within applications to holographic 
superconductors (4l| . A similar problem to that consid¬ 
ered here, for GB theory in vacuum, was studied in (42| 
for AdS spaces and in [43j] for dS spaces. Moreover, ex¬ 
haustive analysis of these space-times in the neutral j44| 
and charged case under Maxwell theory [45| have been 
also performed. The results obtained in this section ex¬ 
tend those for spherical topology horizons k = 1 to the 
general class of admissible NED models supporting finite- 
energy ESS solutions. In our models the equations giving 
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the horizons (l25l) and the extreme black hole radius (l37l) 
receive l\ corrections as 


I/-5 
2 h 


+ r 


D -3 
h 


2 K 2 £ ex {r hl q)\ 

D-2 ) 


D -1 


/ 2 


=(07) 


'■' 2 (rt{r hext r,q) - — --— j = 


while for the extreme black hole mass these corrections 
cancel between themselves and Eg . (13811 remains unmodi¬ 
fied. Concerning the metric function g a {f), its behaviour 
around the center is still dominated by the quantity E(g) 
but the large-r/j behaviour is now given by 


2. l 2 a < 0 

The term under the square root in (f49|) is always posi¬ 
tive and thus l\ is unbounded. The solutions in this case 
are similar as those of subcases (IIV A 2D . (IIVB ID and 
(IIV B 31) . again with a AdS asymptotics r 2 / l^ff Note 
that the typical behaviour of the Reissner-Nordstrom- 
AdS solution is obtained when D > 5,E(g) > 0,Z 2 o < 
Z 2 < 0. It is worth mentioning that in the special case 
Z 2 = 2Z 2 the effective cosmological constant of Ea. fldUl) is 
precisely Z 2 and thus the GB parameter a is identified 
(times a constant) to the effective cosmological constant 
of the theory. 


B. Asymptotically dS solutions (Z 2 < 0) 


g a (r) -A 1 -I- j 2 [ l 



(49) 


and the usual cosmological constant l\ is replaced by the 
effective Z 2 ^^ at large r, defined as 


Z 


2 _ 
eff ~ 


1 - 


Z 


2 

a 




(50) 


However, the signs of l\ and Z 2 ^^ coincide in all cases, 
which leads to the usual correspondence between the sign 
of l\ and the asymptotically (A)dS structures. Moreover, 
the combination of this sign with the one of Z 2 leads to 
four different cases. Let us analyze each case separately, 
using again the sign of E. 


A. Asymptotically AdS solutions ( l\ > 0) 

1. l 2 a > 0 

The term under the square root in (I49[) may become 
negative. Consequently the cosmological constant term 
must be bounded by below as l\ > 2Z 2 for the metric to 
be well defined. When E (q) > 0 the behaviour of the met¬ 
ric for D > 5 becomes similar to subcase (IIV A ID . but 
with an AdS asymptotics r 2 /Z 2 y^. Consequently, when 
D > 5 this space-time is asymptotically similar to the 
Schwarzschild-AdS one. For D = 5 the metric at the 
center is finite and there are again solutions with a single 
nondegenerate horizon if Z 2 < 2E, and two, one (degener¬ 
ate) or zero horizons otherwise, as in the asymptotically 
flat case of subsection IIV Al 

On the other hand, if E(g) < 0 the solutions exhibit a 
BS, becoming similar as those of subcase (IIV B 21) for all 
D , but with an AdS asymptotics as in the previous case. 


This case becomes much more involved due to the in¬ 
terplay between m, e, Z 2 and l\ and, as a consequence, 
there are several additional black hole structures, already 
found in the GR-NED system (see e.g. [I3|)- Let us 
briefly discuss the different possibilities. 


1. l 2 a > 0 

Now l\ is unbounded. If E > 0 and D > 5 the met¬ 
ric diverges to —oo both as r ~ 0 and at r —> oo and 
we find black holes with a extreme (degenerate) hori¬ 
zon with mass m = mh e xtr{q ), naked singularities (m < 
nihextriq)), and two-horizons black holes (m < mhextr ), 
resembling the usual Schwarzschild-dS behaviour. Simi¬ 
lar number of horizons are found when D = 5 if Z 2 < 2E, 
for which the metric at the center Eq. is negative (but 
finite). However, when Z 2 > 2E the metric at the center 
m3 is positive (for A0, Al and A2 fields) and the inter¬ 
play between Z 2 and l\ leads to new kinds of structures 
formed with three types of horizons (inner, outer and cos¬ 
mological) . Indeed now there exist two kinds of extreme 

solutions, with masses m^ xtr (q) < m^ xtr (q) given by 
Eq. (l38l) . In the former case the degenerate horizon is re¬ 
alized through a joining between the inner and outer hori¬ 
zons (with an additional cosmological horizon) represent¬ 
ing the extreme black hole spacetime, while the outer and 
cosmological ones join in the latter, and there is an addi- 
tional inner horizon. For masses such that m > Tn y h ^ xtr (q) 
or m < t (q) there is a solution with a single (cosmo¬ 
logical) horizon. Finally, for m^ xtr (q) < m < m^ xtr (q) 
we are led to black holes with three non-degenerate hori¬ 
zons. This kind of structures are also found in the GB- 
Maxwell case |45| ■ 

On the other hand, if E < 0, as already stated (see 
subsection IIV Bl) the metric shows a BS at r = rs, tak¬ 
ing a (finite) positive (> 1) value there for any D > 4. 
Consequently the structures in terms of horizons are the 
same as in the case D = 5, Z 2 > 2E. 
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2. l 2 a < 0 

Now l\ is bounded as l\ < 21\ and there are three 
subcases for the metric behaviour. When T,(q) > 0 the 
metric at the BS point g a (rs ) = 1 + r| /l^ can be either 
positive or negative, depending on |Z^| < (its deriva¬ 
tive there is always positive): the former case leads to a 
single-horizon black hole while in the latter there can be 
black holes with two horizons (event and cosmological), 
a single one (extremal, the event and the cosmological 
horizons join) or naked singularities, depending on the 
mass according to m = mhextr(q )■ Note that, as in the 
previous cases, l\ and r| are not independent, being re¬ 
lated through equation (l27l) . 

When £ (q) < 0 and l^ 0 < /^ < 0, as already 
seen, the metric diverges around the center when D > 
5 while for D = 5 it takes a positive finite value 

g a (0) —>• 1 — -— — > 1, but the structure of hori¬ 
zons in both cases is similar as that found in subcase 
(IV B ID with £(g) < 0: extreme solutions m = m^ xtr 
(degenerate inner+outer horizon and cosmological hori- 

zon) and m = frv h ^ xi . r (inner horizon and degenerate 
outer+cosmological horizon), black holes with three hori¬ 
zons m^ xtr (q) < m < m^ xtr (q) and solutions with 
a single cosmological horizon for m < ni^ xtr (q) or 

to > ^hlxtril)- Finally the case £(g) < 0, 1% < Z„ 0 < 0 
shows two branch singularities and, since the region be¬ 
yond the outer BS radius is the physically relevant, the 
structure of horizons is similar as in case £ < 0, 1^ > 0. 


VI. CONCLUSIONS 

In this paper the families of NED models, constrained 
by several physical admissibility requirements, and sup¬ 
porting ESS solutions whose energy in flat space is finite, 
were considered in the framework of Gauss-Bonnet the¬ 
ory. This is the simplest nontrivial extension of General 
Relativity incorporating higher-order curvature invari¬ 
ants and leading to second-order field equations. The fi¬ 
nite character of the energy of the ESS field is established 
according to the asymptotic and central-field behaviour 
of the ESS field in absence of gravity. In the latter case 
two held behaviours compatible with the finite-energy re¬ 
quirement were determined: one divergent at the origin 
and the other one attaining a finite value there. 

With these results we have shown that the associated 
gravitating structures for these NED models, when cou¬ 
pled to the Einstein-Gauss-Bonnet action, can be qualita¬ 
tively characterized in terms of the relation between the 
mass parameter m and the ESS energy e(q). These struc¬ 
tures fall into two classes according to the Gauss-Bonnet 
parameter l ^, one defining solutions whose domain of ex¬ 
istence is the whole spacetime, and another for which the 
solutions are not defined everywhere, showing a singular¬ 
ity at a finite radius. In the former case the nature of 


the ESS field at the center (A0, Al, A2) and the space- 
time dimension D critically affects the metric behaviour. 
Indeed, the case D = 5 was shown to posses a different 
structure as compared with the remaining cases D > 5. 
In addition, when a cosmological constant term is intro¬ 
duced in the system, then the relation between the mod¬ 
ule and sign of l\ and l\ determines the effective value of 
the cosmological constant of the corresponding asymp¬ 
totically (Anti-)de Sitter space. While in the former case 
the number and kind of gravitating structures remain 
unmodified, in the latter new structures may appear, in¬ 
cluding black holes with three non-degenerate horizons, 
or black holes with a both a degenerated and a nondegen¬ 
erate horizon. The analysis performed here goes beyond 
previous results obtained in the literature corresponding 
to several particular Lagrangian densities [HJ, and ex¬ 
tends them to the family of physically admissible grav¬ 
itating NED models with finite-energy ESS solutions in 
the context of Gauss-Bonnet theory. 

Although we have only analyzed here black holes with 
the usual spherical topology for the horizon, it is possible 
to extend these results to the cases of negative constant 
or zero curvature hypersurfaces. Also a suitable exten¬ 
sion of this work would be the analysis of the thermody¬ 
namic features of these solutions, concerning also those 
with the different horizon topologies mentioned above. 
Indeed, it is well known that while usually entropy of 
black holes equals one quarter of the horizon radius, this 
property does not hold, in general, for higher-order cur¬ 
vature theories [46}- Thus, it would be worth studying 
such features, in a similar way as is done in [42], j43| for 
(A)dS-GB theory in vacuum and with a Maxwell field. 
It is expected that, analogously to the geometric study 
carried out in this paper, thermodynamics of the corre¬ 
sponding solutions will also depend on a few data and, 
therefore, a similar systematic analysis will also be possi¬ 
ble. Such investigations, for non-vanishing cosmological 
constant term, could be of potential relevance within the 
framework of gauge/gravity dualities. At this point of 
our research, and in absence of specific settings on the 
CFT side, it is difficult to foresee what kind of particular 
NED theory would be useful on the gravity side. This 
further supports the interest on methods like those de¬ 
veloped here to characterize the generic geometric and 
thermodynamic behaviour of physically consistent NED 
theories in asymptotically AdS space-times, which could 
prove useful in the future to match particular CFT set¬ 
tings. 

To conclude, since Gauss-Bonnet and, more generally, 
Lovelock gravities, are at the crossroad of the metric and 
Palatini formulations Q, the study of these higher-order 
gravity theories and their solutions may reveal useful in¬ 
formation on the effective approaches to quantum grav¬ 
ity. A final remark concerns the fact that GB theories 
coupled to admissible NED models are unable to pro¬ 
vide a solution free of curvature singularities everywhere, 
a reminiscent situation as that of GR, where the cen¬ 
tral singularity is only avoided (for purely electric fields) 
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through adhoc unphysical choices of the NED Lagrangian 
density 0 - This motivates further investigations on the 
avoidance of curvature singularities in other scenarios of 
modified gravity. 
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